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Abstract: 

The radiative corrections of the strong and electroweak interactions are calculated for 
the Higgs-boson decays H — > WW/ZZ — > 4/ with semileptonic or hadronic four-fermion 
final states in next-to-leading order. This calculation is improved by higher-order cor- 
rections originating from heavy-Higgs-boson effects and photonic final-state radiation off 
charged leptons. The W- and Z-boson resonances are treated within the complex-mass 
scheme, i.e. without any resonance expansion or on-shell approximation. The calculation 
essentially follows our previous study of purely leptonic final states. The electroweak cor- 
rections are similar for all four-fermion final states; for integrated quantities they amount 
to some per cent and increase with growing Higgs-boson mass Mr, reaching 7-8% at 
Mh ~ 500 GeV. For distributions, the corrections are somewhat larger and, in general, 
distort the shapes. Among the QCD corrections, which include corrections to interference 
contributions of the Born diagrams, only the corrections to the squared Born diagrams 
turn out to be relevant. These contributions can be attributed to the gauge-boson de- 
cays, i.e. they approximately amount to a s /ir for semileptonic final states and 2a s /7f for 
hadronic final states. The discussed corrections have been implemented in the Monte 
Carlo event generator Prophecy4f|3 
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1 Introduction 



The startup of the Large Hadron Collider (LHC) in 2007 will open up a new era in 
particle physics. One of the main tasks of the LHC will be the detection and the study of 
the Higgs boson. If it is heavier than 140 GeV and behaves as predicted by the Standard 
Model (SM), it decays predominantly into gauge-boson pairs and subsequently into four 
light fermions. From a Higgs-boson mass M H of about 130 GeV up to the Z-boson-pair 
threshold 2Mz, the decay signature H — > WW* — > 21eptons + missing pp [1] has the 
highest discovery potential for the Higgs boson at the LHC [2]. For higher Higgs-boson 
masses, the leading role is taken over by the "gold-plated" channel H — > TIL — > 4 leptons, 
which will allow for the most accurate measurement of M H above 130 GeV [3] . More 
details and recent developments concerning Higgs-boson studies at the LHC can be found 
in the literature [4,5]. At a future e + e~ linear collider [6], the decays H — > 4/ will enable 
measurements of the H — > WW/ZZ branching ratios at the level of a few to 10% [7]. 

At the LHC, owing to the huge background of strongly interacting particles, the most 
important decay modes in H — > WW/ZZ — > 4/ are those with leptons in the final state. 
Therefore, most analyses are based on them. However, also final states involving quarks 
can be useful owing to their larger branching fractions. For decays involving intermediate 
W bosons these provide better kinematical information since they involve less neutrinos. 
For instance, it has been found that in the vector-boson-fusion channel the decays H — > 
WW — > l^ujj can provide complimentary evidence in the intermediate Higgs-mass range 
140 GeV < M H < 200 GeV [2,8,9] and constitute a good potential discovery channel in the 
medium-high Higgs-mass range Mh > 300 GeV [4] . At a linear collider, the hadronic and 
semileptonic final states are even more important since they allow for a full reconstruction 
of the Higgs-boson decay H — > WW [10]. 

A kinematical reconstruction of the Higgs-boson decays H — > WW — > 4/ and the 
suppression of the corresponding backgrounds requires the study of distributions and the 
use of cuts defined from the kinematics of the decay fermions. In addition, the verification 
of the spin and of the CP properties of the Higgs boson relies on the study of angular, 
energy, and invariant-mass distributions [11]. These tasks require a Monte Carlo generator 
for H — > WW/ZZ — > 4f. Since the effects of radiative corrections, in particular real- 
photon or gluon radiation, are important, a Monte Carlo generator including all relevant 
corrections is needed. 

The progress in the theoretical description of the decays of a SM Higgs boson into W- or 
Z-boson pairs has, for instance, been summarized in Ref. [12]. Until recently, calculations 
for off-shell vector bosons were only available in lowest order [13], and radiative corrections 
were known only in narrow- width approximation (NWA) [14], i.e. for on-shell W and 
Z bosons. In this case, also leading two-loop corrections enhanced by powers of the top- 
quark mass [15] or of the Higgs-boson mass [16] have been calculated. However, near 
and below the gauge-boson-pair thresholds the NWA is not applicable, so that only the 
lowest-order results existed in this Mh range. 

In a recent paper [12] we have presented results for the complete electroweak (EW) 
0(a) corrections including some higher-order improvements to the Higgs-boson decays 
H — > WW/ZZ — > 4 leptons. First results of this calculation had already been presented 
at the RADCOR05 conference [17]. At this conference also progress on an independent 
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calculation of the electromagnetic corrections to H — > TIL — > 4 leptons has been reported 
by Carloni Calame et al. [18]. The analytic results demonstrated in Ref. [12] are also 
valid for quarks in the final state. In this paper we supplement this calculation by the 
corresponding QCD corrections. We introduce a classification of the QCD corrections and 
describe their calculation. The QCD corrections have been implemented into the Monte 
Carlo generator Prophecy4f, and numerical results have been produced. These include 
the partial widths for various semileptonic and hadronic channels as well as different 
invariant-mass and angular distributions for semileptonic final states. 

The paper is organized as follows: In Section [2]we describe the setup of our calculation. 
Section [3] contains a classification of the QCD corrections and provides analytic results 
for the virtual and real QCD corrections. Numerical results are presented in Section HI 
and our conclusions are given in Section [5j 

2 Setup of the calculation 

We consider the processes 

H(p) — fiih,^) + f 2 (k 2 ,a 2 ) + f 3 (k 3 ,a 3 ) + Uih,^) + [j/g(k,X)}, (2.1) 

where fa stands for any lepton, / = e, //, r, u e , u^, u T , or for any quark of the first two 
generations, q = d, u, s, c. We do not include final states with bottom or top quarks. 
The momenta and helicities of the external particles are indicated in parentheses. The 
helicities take the values <Ji = ±1/2, but we often use only the sign to indicate the helicity. 
The masses of the external fermions are neglected whenever possible, i.e. everywhere 
but in the mass-singular logarithms. We always sum over the four light quarks of the 
first two generations in the final state and set the CKM matrix to the unit matrix. 
This approximation ignores quark mixing with the third generation, which is, however, 
negligible. 

The calculation of the EW corrections has already been described in Ref. [12], where 
results for purely leptonic final states have been discussed. Here we briefly repeat the 
salient features of the evaluation of virtual one-loop and real-photonic corrections. 

The calculation of the one-loop diagrams has been performed in the conventional 
't Hooft-Feynman gauge and in the background-field formalism using the conventions of 
Refs. [19] and [20], respectively. 

For the implementation of the finite widths of the gauge bosons we use the complex- 
mass scheme, which was introduced in Ref. [21] for lowest-order calculations and gener- 
alized to the one-loop level in Ref. [22]. In this approach the W- and Z-boson masses are 
consistently considered as complex quantities, defined as the locations of the propagator 
poles in the complex plane. The scheme fully respects all relations that follow from gauge 
invariance. A brief description of this scheme can also be found in Ref. [23]. 

The amplitudes have been generated with FeynArts, using the two independent 
versions 1 and 3, as described in Refs. [24] and [25], respectively. The algebraic evaluation 
has been performed in two completely independent ways. One calculation is based on an 
in-house program written in Mathematica, the other has been completed with the help of 
FormCalc [26]. The amplitudes are expressed in terms of standard matrix elements and 
coefficients, which contain the tensor integrals, as described in the appendix of Ref. [27]. 
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The tensor integrals are evaluated as in the calculation of the corrections to e + e~ — > 4/ 
[22,28]. They are recursively reduced to master integrals at the numerical level. The 
scalar master integrals are evaluated for complex masses using the methods and results 
of Ref. [29]. UV divergences are regulated dimensionally and IR divergences with an 
infinitesimal photon mass. Tensor and scalar 5-point functions are directly expressed in 
terms of 4-point integrals [30]. Tensor 4-point and 3-point integrals are reduced to scalar 
integrals with the Passarino-Veltman algorithm [31] as long as no small Gram determinant 
appears in the reduction. If small Gram determinants occur, two alternative schemes are 
applied [32]. One method makes use of expansions of the tensor coefficients about the 
limit of vanishing Gram determinants and possibly other kinematical determinants. In 
the second, alternative method we evaluate a specific tensor coefficient, the integrand 
of which is logarithmic in Feynman parametrization, by numerical integration. Then 
the remaining coefficients as well as the standard scalar integral are algebraically derived 
from this coefficient. The results of the two different codes, based on the different methods 
described above are in good numerical agreement. 

Since corrections due to the self-interaction of the Higgs boson become important 
for large Higgs-boson masses, we have included the dominant two-loop corrections to the 
decay H — > VV proportional to G^M^ in the large-Higgs-mass limit which were calculated 
in Ref. [16]. 

The matrix elements for the real-photonic corrections are evaluated using the Weyl- 
van der Waerden spinor technique as formulated in Ref. [33] and have been checked against 
results obtained with Madgraph [34]. The soft and collinear singularities are treated 
both in the dipole subtraction method following Ref. [35] and in the phase-space slicing 
method following Ref. [36]. For the calculation of non-collinear-safe observables we use 
the extension of the subtraction method introduced in Ref. [37]. Final-state radiation 
beyond 0(a) is included at the leading-logarithmic level using the structure functions 
given in Ref. [38] (see also references therein). 

The phase-space integration is performed with Monte Carlo techniques. Prophecy4f 
employs a multi-channel Monte Carlo generator [39] similar to the one implemented in 
RACOON WW [21] and C0FFER77 [37,40]. Our second code uses the adaptive integration 
program VEGAS [41]. 

3 QCD corrections for H — > 2q2l and H — > 4q 
3.1 Classification 

A proper classification of QCD corrections is achieved upon considering possible contri- 
butions to the squared lowest-order amplitude. The amplitude itself receives contributions 
from one of the two tree diagrams shown in Figure [T] or from both. Thus, the square of 
this amplitude receives contributions from cut diagrams of the types depicted in Figure [21 
Type (A) corresponds to the squares of each of the Born diagrams, type (B) to their 
interference if two Born diagrams exist. 

After this preliminary consideration we define four different categories of QCD correc- 
tions. Examples of cut diagrams belonging to these categories are shown in Figure El the 
corresponding virtual QCD correction diagrams are depicted in Figure HI 
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Figure 1: Possible lowest-order diagrams for H — > 4f where V, V = W, Z. 




Figure 2: Types of cut diagrams contributing in lowest order. 



(a) QCD corrections to gauge-boson decays comprise all cut diagrams resulting from 
diagram (A) of Figure [2] by adding one additional gluon. Cut diagrams in which 
the gluon does not cross the cut correspond to virtual one-loop corrections, the one 
where the gluon crosses the cut correspond to real-gluon radiation. Note that cut 
diagrams in which the gluon connects the two closed quark lines identically vanish, 
because their colour structure is proportional to Tr(A /l )Tr(A ,! ') = 0, where X h is 
a Gell-Mann matrix. Thus, the only relevant one-loop diagrams in this category 
are gluonic vertex corrections to a weak-boson decay, as illustrated in the first 
diagram of Figure HI the real corrections are induced by the corresponding gluon 
bremsstrahlung diagrams. 

If a weak-boson decay is fully integrated over its decay angles, the resulting QCD 
correction of the considered type simply reduces to the well-known factor a s /ir for 
a hadronically decaying vector boson. 

(b) QCD corrections to interferences comprise all cut diagrams resulting from diagram 
(B) of Figure [2]by adding one additional gluon, analogously to the previous category. 
Relevant one-loop diagrams are, thus, vertex corrections or pentagon diagrams, as 
illustrated in the first two diagrams of Figure HI 

(c) Corrections from intermediate qqg* states are induced by loop diagrams exemplified 
by the third graph in Figure HI The remaining graphs are obtained by shifting the 
gluon to different positions at the same quark line and by interchanging the role of 
the two quark lines. Thus, the diagrams comprise not only box diagrams but also 
vertex diagrams. They do not interfere with Born diagrams with the same fermion- 
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Figure 3: Categories of cut diagrams contributing to the QCD corrections. 




Figure 4: Basic diagrams contributing to the virtual QCD corrections for H — > 4/ where 
V = W, Z and q — d, u, s, c, b, t. The categories of QCD corrections, (a)-(d), to which the 
diagrams contribute are indicated. 
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number flow because of the colour structure, i.e. in 0(a B ) they only contribute if 
two Born diagrams exist. 

Owing to the intermediate qqg* states, the squared diagrams of this category actually 
correspond to (collinear-singular) real NLO QCD corrections to the loop-induced 
decay H — ► qqg, where q is a massless quark. Here we consider only the interference 
contributions of the loop diagrams of this category with the lowest-order diagrams 
for the decay H — ► VV — > 4q, resulting in a UV and IR (soft and collinear) finite 
correction. 

(d) Corrections from intermediate g*g* states are induced by diagrams exemplified by 
the fourth graph in Figure 0J There are precisely two graphs with opposite fermion- 
number flow in the loop. Again, owing to the colour structure (see also below), 
these diagrams do not interfere with Born diagrams with the same fermion-number 
flow, i.e. the existence of two Born diagrams is needed. 

Owing to the intermediate g*g* states, the squared diagrams of this category actually 
correspond to (collinear-singular) real NNLO QCD corrections to the loop-induced 
decay H — > gg. The considered interference contributions of the loop diagrams of 
this category with the lowest-order diagrams for the decay H — > VV — *> 4q, however, 
again yield a UV and IR (soft and collinear) finite correction. 

From the classification, it is clear that category (a) exists for all final states involving 
quarks, while categories (b), (c), and (d) are only relevant for the hadronic decays H — > 
qqqq and H — > qqq'q', where q and q' are weak-isospin partners. Categories (a), (b), and 

(c) give rise to contributions to the decay widths that are proportional to a 3 a s , while type 

(d) yields a contribution proportional to a 2 a 2 . 

We do not consider the process H — > 4 jets in general but only the contributions 
via virtual EW gauge-boson pairs, i.e. we assume that the gauge-boson resonances are 
isolated by experimental cuts. For the more inclusive decay H — > 4 jets, also diagrams 
without intermediate EW gauge bosons, where the Higgs boson couples to gluons via 
heavy-quark loops, become important. Using an effective Hgg coupling, the calculation 
of the corresponding QCD one- loop matrix elements has been described in Ref. [42], but 
the full NLO QCD prediction for H — > 4 jets including these effects is not yet available. 

3.2 Virtual corrections 

In the evaluation of the one-loop QCD diagrams, which are illustrated in Figure HI 
the fermion spinor chains are separated from the rest of the amplitude by introducing 
52 standard matrix elements jCi® cd ' a7 , as defined in Eq. (3.2) of Ref. [12], where the 
indices a and r indicate the chiralities in the spinor chains of the fermion pairs f a fb and 
fcfdi respectively. Furthermore, the colour structure is extracted by defining the colour 
operators 

Cf cd = 5 CCh ® 5 C Crl , Cf cd = -^—\ h c c ®\ h cc = —\ h c c ®\ h cc (3.1) 

with the Gell-Mann matrices X h , the colour index h of the gluon, and the colour indices 
c a ,b,c,d of the quarks. For external leptons the corresponding colour index trivially takes 
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only one value, and the operator C 2 , of course, appears only for four-quark final states. 
Using this notation, the generic lowest-order amplitude in colour space reads 

•^o,cadccCd d (k a , kt,, k c , ka) = C 1 J\4 b d (k a , kb, k c , kd), (3-2) 

where j\4^ v ^ aCTbCTcCTd j s the colour-stripped generic lowest-order amplitude defined in 
Eq. (2.7) of Ref. [12]. Obviously, this notation generalizes to the generic EW one-loop 
amplitudes (i.e. without gluon exchange) introduced in Eq. (3.3) of Ref. [12] 



-^EW ,c a c b c c c d — b l -^'EW 



abed \ AVV^aObVcOd 



13 



G l Z^M 



abcd,a a cTc iCAabcd,a a a c x x 
EW,i jvl i °cr a ,-a b Oa c ,-a d , 



(3.3) 



where j£\°: c >° a(Tc denote the standard matrix elements and F^ff aUc are Lorentz-invariant 
coefficient functions. In the generic amplitudes the superscript U VV" indicates the com- 
mon fermion-number flow, which corresponds to the decays V — > f a fb and V — > f c fd- The 
one-loop QCD amplitude, which involves gluon exchange, receives contributions from both 
colour operators; in colour space we define 



j^/V,<T a <Jb&c<Jd \^ ^yabed \ yfVV,a a <JbCr c <jd 



Efiabcd ka V 
13 

j^VV,a a cr b a c ad jpabcd,a a cr c j^abcd,a a cr c ^ 



V F a 

i=l 



-<?l .o- c ,—a d i 



(3.4) 



where the M^ci) a J bUcUd are colour-stripped amplitudes. 

From the generic matrix elements AYY'Z a ° b ? cCrd (n = 0, 1) the matrix elements A° a ° b ° c ?i 
for the specific processes are constructed as in Eqs. (2.11)-(2.14) of Ref. [12]. The index 
n = 1 collectively represents the sum EW + QCD of EW and QCD one- loop contributions. 
We denote different fermions by / and F, and their weak-isospin partners by /' and F' 
(/ 7^ F,F'). For purely hadronic final states the quarks are denoted by q and their 
weak-isospin partners by q' . Thus, we obtain: 



H - ffFF: 



(3.5) 



H -> ff'FF': 



^nx7c2C3C 4 (^l' ^2, h, ki) 



A 



WW,<ti 0-20-30-4 

n,ClC2C3C4 



(fa, k 2 , k 3 , k 



(3.6) 



H 



qqqq: 



An]ciC2 3 c°c 4 (ki, k 2 , k 3 , k 4 ) — A^^"J i i U4 (k l) k 2 , k 3 , k 4 ) 

aZZ, 01 040302 (r, u u u \ 
"^n,c\CACzC<> 1^1; ft 4, ^3) «-2j, 



(3.7) 



1 We note that the generic colour-stripped EW one-loop amplitude M^ f f aC ' h ' J '' ad was denoted 



abed, err 



in Eq. (3.3) of Ref. [12]. 
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• H — > qqq'q': 

A a n°iC2c z c4,{ki-, k 2 , k 3 , k 4 ) = J^c°c2?°ct Ci (ki, k 2 , k 3 , k 4 ) 

-^:z 4 r 2 (h,h,k 3} k 2 ). (3.8) 

The relative signs between contributions of the basic subamplitudes to the full matrix 
elements account for the sign changes resulting from interchanging external fermion lines. 

Since the lowest-order amplitudes only involve the colour operators C| 234 and C 1432 , 
the following colour sums appear in the calculation of squared lowest-order amplitudes 
and of interferences between one-loop and lowest-order matrix elements: 

X[ A) = ^(Cf^Cf™ 1 ) = N c fa N c fc , x\ A) = Y,(Cf cd *C% bcd ) = 0, 

{Ci} {Ci} 

X[ B) = ^{C^Cf^) = N c fa , xi B) = Y^iC^Cf*) = N c fa , (3.9) 

{Ci} {Ci} 

where J2{ Ci } stands for the sum over the colour indices c a , Cf,, c c , Cd, and Nj is the colour 
factor for a fermion /, which is 1 for leptons and 3 for quarks. 

Squared Born diagrams, as illustrated in type (A) of Figure [U are proportional to x[ A \ 
lowest-order interference diagrams of type (B) are proportional to x[ . The situation 
is analogous for all one-loop diagrams without gluons. By definition, category (a) of 
the gluonic diagrams comprises all one-loop QCD corrections proportional to x[ . In 
category (b), the vertex corrections are proportional to x[ B ^ and the pentagons to X% . 
Categories (c) and (d) receive only contributions from X 2 ; interferences of one-loop 
diagrams like (c) and (d) in Figure H] with Born diagrams of the same fermion-number 
flow vanish because of X 2 A ^ = 0. 

Finally, we obtain the following for the one-loop corrections to the squared matrix 
elements: 

• H -> ffFF: 

E 2 Re {Al%%lX A%™™ } = 2 Re {N c f N c F M f>-™ *A^-™ } , ( 3 . 10) 

{Ci} 

• H -> ffFF': 

E9 Ro / AO\01 °"3 °"4 * AmCT2 T3CT4 \ _ 9 Tj / Arc Arc « >WW,<Tl(T2Cr30-4 * « A WW,<TlCT21T3e7-4 \ /O II - ! 
Zn - e \- /l 0,ciC2C 3 C4 • /l l,ciC2C 3 C4/ — Zn - e \ JV / JV F- /w ■ /W EW+QCD(a) /> \°- lL ) 

{Ci} 

• H — > qqqq: 

EOTDpf ^0-10-20-30-4* ^JCTlO-203CT4 \ 
^ 1VC \ ,/ n),ClC2C3C4 •^ 1 l,ClC2C3C4j 

id} 

o T>„( A / |ZZ,oi02CT304 * \( A7"C\2 A^ZZ,0-lO-203CT4 t»tc I /ZZ,CTl04CT302 1 \ 

— Zrte^yviQ |_l iV aJ ■ /vt EW+QCD(a) JV <? ■ /w EW+QCD(b)+QCD(c)+QCD(d)J J 

+ (g(^ 2 ,cr 2 ) ^g(fc 4 ,cT 4 )), (3.12) 
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H — > qqq'q': 

E9 Rp / J"l ff 2ff3"4 * 40-10-20-30-4 \ 
Z IXtJ ^0,ClC 2 C 3 C4 • /i l,CiC2C3C 4 / 

opj A / |ZZ,CTiCT 2 03Cr4 * [/ A7-CN2 «/ZZ, 01 020304 »rc »^ WW, 01 040302 

— Ztte| yvi [\ ly q) • /W EW+QCD(a) JV g - /W EW+QCD(b) 



. ;WW, 01O4O3O2 * / l\TC\2 « / WW, 01040; 

"T" JVl \ Iy q) • /W EW+QCD(; 

— N C A / | ZZ . <J l <T 20-304 1 1 /o 10\ 

iv g - /w EW+QCD(b)+QCD(c)+QCD(d)J /• 1,-J.iO; 

Due to the electric charge flow, categories (c) and (d) only exist if there are corresponding 
diagrams with intermediate Z bosons. That is why there are no terms _/\ / ^ zz ' <T10 " 20 "3o-4 * x 
^QOT(c)+QCD(d) ■ Note that in the notation we have suppressed the momentum arguments 
which, however, can be trivially restored, because the permutation of momenta ki is the 
same as for the polarizations Oi in each amplitude. 

3.3 Matrix element for real-gluon emission H — > 4fg 

The real-gluonic corrections are induced by the process 

H(p) — fi(h, a,) + f 2 (k 2 , <r 2 ) + Hh, a 3 ) + f 4 (h, a 4 ) + g(k, A), (3.14) 

where the momenta and helicities of the external particles are indicated in parentheses. 

The matrix elements for this process can be constructed from the matrix elements for 
the photon radiation process A4^ aabacadX (Q a , Qb, Q c , Qd, k a , kb, k c , k d , k), which have been 
explicitly given in Ref. [12]. Here, Q a ,b,c,d denote the electric charges of the fermions. The 
generic amplitudes read 



<^ 2 A ^WW, 01040302 



AVV,cr a a b a c a d \,h/i i i i i\ 
*H-R,c a c h c c c d \K a , Kb, Kc, Kd, K) 

e I 2 



g,C a C b C c C d 

1 ' 1 ^icA^ M V' aaabacadX ^ hO,0,k a ,h,k c ,k d , k) 



+ ^t/ Cl ^/X^ WdA ( ' 0» !» !» **» k, ^» fc )}> ( 3 - 15 ) 

where g s is the strong coupling constant, and V = Z, W for Z-mediated and W-mediated 
decays, respectively. The symbols 5f iq are equal to one if /j is a quark and zero otherwise. 

From the generic matrix element Ag^^ c b ^ d cadX,h (k a , k b , k c , k d) k) the matrix elements for 
the specific processes can be constructed as follows. As above, we denote different fermions 
(/ 7^ F, F') by / and F, and their weak-isospin partners by /' and F', respectively: 

• H - ffFFg: 

Kc^\ k ^ ^ fc 3, h, k) = A^Zr^ih, k 2 , k 3 , h, k), (3.16) 

• H - ff'FF'g: 

AZ?c? c Ict h (h, k 2 , fcs, h, k) = Af^Z^'Hh, k 2 , k 3 , h, k), (3.17) 



9 



H -> qqqqg: 



^Z7c? c ^ x ' h (h, h, h, fc 4 , fc) = A*%lZr x ' h (h, k 2 , k 3 , fc*, fc) 

- ^\zzr x ' h (^ fcs. **, (3-i8) 

• H — > qqq'q'g: 

^^c x ' h (h, h, fc s , fc 4 , fc) = ^s«r 4A,fc (*i. **> fc 3 , fc 4 , fc) 

- CcSr^d. **, fc 2 , fc). (3.19) 

The relative signs between contributions of the basic subamplitudes to the full matrix 
elements account for the sign changes resulting from interchanging external fermion lines. 
Squaring the amplitudes and summing over the colour degrees of freedom, we have 



H - ffFFg: 

E \^Z^ h (h,h,h,h,k) 

{ci},h 



2 = ^N c f N c F ^ 
3 f F a 



x 



>/5 



Mf^ a ^\l, 1, 0, 0, h, fc 2 , fc 3 , fc 4 , fc) 



+ 5 Fq M^'^^iO, 0,1,1, fci, fc 2 , fc 3 , fc 4 , fc) 



(3.20) 



• H - ff'FF'g: 



E I^Ss^'^fci, fc 2 , fcs, fc 4 , fc) = ^n f ^ 

{c;},/l 



2 4 



X 



M^' 4 ^ 1 ' °> °> fci> fc2 ? fc s, fc 4 , fc) 



+ 5 Fq M^ w ' aia2a3(T4X (0, 0, 1, 1, fc 1; fc 2 , fc 3 , fc 4 , k) 



(3.21) 



• H -> qqqqg: 



E |^c? C Xc 4 A ' ft (fcl,fc2,fcs,fc 4 ,fc) = gTO 2 " 
{ci},h 



2 4 



a 



x 



+ 
+ 
+ 



-Mf ' CTlCT2a3<J4A (l, 1, 0, 0, fc 1; fc 2 , k 3 , k 4 , k)\ 
Mf^ a ^\0, 0, 1, 1, fci, fc 2 , fc 3 , fc 4 , fc) 
M; z ' W4W2A (l, 1, 0, 0, fc x , fc 4 , fc 3 , fc 2 , fc) 
^f 51W3a2A (0, 0, 1, 1, h, k 4 , fc 3 , fc 2 , fc) 

r (Mf^^\l, 1, 0, 0, h, k 2 , fc 3 , fc 4 , fc) 
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+ Mf' aia2asa4X (0, 0,1,1, fa, fa, fa, fa, k 
x (M^w^il, 1, 0, 0, fa, fa, fa, fa, k) 

+ Mf' aia4asa2X (0, 0,1,1, fa, fa, fa, fa, k 



• H 



qqq'q'g: 



£ 

{a},h 



\^z? c ?cZ x ' h (h, fa, k 3 , fa, k) 2 = t(K) 2 - 



a 



x 



(3.22) 



Mf' aia2a3a4X (l, 1, 0, 0, fa, fa, fa, fa, k)\ 

+ |tw^ z ' CTiCT2CT3CT4A (0, 0, 1, 1, fa, fa, fa, fa, k)\ 2 
+ \M^ W ^^ 2 \1, 1, 0, 0, fa, fa, fa, fa, k) 
+ M^ w ^^ 2 \0, 0, 1, 1, fa, fa, fa, fa, k) 

r (Mf^^\l, 1,0,0, fa, fa, fa, fa, k) 

+ Mf^ c ^\0, 0,1,1, fa, fa, fa, fa, k))* 
x (^7' ffIffWff2A (l, 1, 0, 0, fa, fa, fa, fa, k) 

+ M™'" 1 ' 4 ^ 2 ^, 0, 1, 1, fa, fa, fa, fa, fc))| . (3.23) 
The contribution T g of the radiative decay to the total decay width is given by 



WH J { Cl },h K},a=±i 



2Mi 



where the phase-space integral is defined by 



C2C3C4 I 



(3.24) 



(3.25) 



4 Numerical results 

4.1 Setup and input 

We use the scheme, i.e. we define the electromagnetic coupling by q; Gm = 
V2G^M^,(1 — Myj/MD/n. Our lowest-order results include the C(a)-corrected width 
of the gauge bosons. In the QCD corrections we uniformly take a fixed value for 
ai s = a s (Mz) = 0.1187 everywhere, because the only numerically relevant part (see below) 
of the QCD correction is the one connected with the hadronic decay of a W or a Z boson, 
where the scale is fixed by the intermediate gauge-boson decay. More details about the 
setup and all input parameters are provided in Ref. [12]. 
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In our approach the final states involve either four fermions (from lowest order and 
virtual corrections), four fermions and a photon (from real-photonic corrections), four 
fermions and a gluon (from real-gluonic corrections), or four fermions and one or more 
photons collinear to an outgoing lepton (from the structure functions describing multi- 
photon final-state radiation). In particular there are no events containing both photons 
and gluons. Moreover, we only consider semileptonic final states in the distributions. For 
these distributions, a photon and gluon recombination is performed as follows. In events 
with a real photon, as in Ref. [12] the photon is recombined with the (in this sense) 
nearest charged fermion if the invariant mass of the photon-fermion pair is below 5 GeV. 
This, in particular, implies that all photons collinear to a lepton are recombined with the 
corresponding lepton if the recombination is switched on (as always done in the results 
of this paper), i.e. the higher-order effects from photonic final-state radiation described 
in Section 4.3 of Ref. [12] fully cancel out in this case. In the case of real-gluon radiation 
we force a 2-jet event. This is achieved by always recombining the two partons of the 
qqg system that yield the smallest invariant mass. Invariant masses and angles are then 
defined by the 4-momenta of the recombined pair and the remaining partons. 

We always sum over the quarks of the first two generations, q — u, d, c, s, and over 
the three neutrinos in the final states and we consider the final states eegg, vvqq, evqq, 
and qqqq. Since we consistently neglect the masses of external fermions and average over 
polarizations, we can express the partial widths of these final states as 

rH->ccgg = 2rH-+ e -e+uu + 2r H ^ c - c + dd , 

^U^uvqq — 6rH^j, c p e uu + ^B_^u c D c dAi 
^Yi^cvqq = 4F H ^ e -p cUd , 

^H^qqqq = Th^iiAcc + TH^ddss + 2Fh^uuss + 2F H ^ udsE 

+ 2r H ^ uddfl + 2r H ^ U fluu + 2r H _ >dddd . (4.1) 

Note that Tn^ euqq includes both electrons and positrons in the final state. The partial 
widths with muons in the final state can be classified in the same way and are equal to 
those with the muons replaced by electrons, because no dependence on the final-state 
fermion masses remains for these inclusive quantities. 

The results for the partial decay widths in the plots are calculated using 10 7 Monte 
Carlo events, while all other results (decay widths in the table and distribution plots) are 
obtained with 5 x 10 7 events. In the presented results, soft and collinear divergences are 
treated with the dipole-subtraction method and have been checked by applying the phase- 
space slicing method. For the latter method more Monte Carlo events are needed for an 
accuracy at the per-mille level, because the energy and angular cuts in this method have 
to be chosen small enough rendering the real corrections and the analytically integrated 
soft and collinear singular contribution (which compensate each other) very large. In 
both methods it is possible to evaluate the virtual corrections (rendered finite by adding 
the soft and collinear singularities from the real corrections) less often than the lowest- 
order matrix elements, because the virtual corrections and also their statistical error 
are smaller. We evaluate the EW virtual corrections only every 100th time and the 
virtual QCD corrections only every 20th time. This procedure reduces the run-time of 
the program while maintaining the size of the overall statistical error. 
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4.8 




OCT) 
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P,W 

i-J vv 


1.8752(4) 
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550.47(7) 
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QCD 


1.9511(7) 


7.5 


159.83(4) 
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566.2(1) 


7.6 




LO 


1.8150(4) 




148.59(2) 
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Table 1: Partial decay widths Th^4/ in lowest order (LO), including 0(a) and 0{G'j 1 M^) 
EW corrections, 0(a s ) QCD corrections, and the sum of EW and QCD corrections (cor- 
rected) and corresponding relative corrections S for semileptonic and hadronic decay chan- 
nels and different Higgs-boson masses. 

4.2 Results for partial decay widths 

In Table [1] we show the partial decay widths of the Higgs boson for semileptonic and 
hadronic final states for different values of the Higgs-boson mass. We list the lowest-order 
(LO) predictions and the predictions including the complete EW 0(a) plus 0{G 2 ^M^) 
corrections and the 0(a s ) QCD corrections. In addition we give the predictions including 
only the EW corrections and only the QCD corrections. In all cases we provide also the 
relative corrections 5 = T/T — 1 in per cent. The statistical errors of the phase-space 
integration are given in parentheses. The size of the EW corrections is very similar to 
the size of the corresponding corrections for leptonic final states discussed in Ref. [12]. 
Since the QCD corrections mainly arise from vertex corrections and since we consider the 
integrated partial widths, the QCD contribution roughly amounts to a s /ii for semileptonic 
final states and 2a s /ii for the hadronic final state. The sum of EW and QCD corrections 
thus rises to 5-14%. 
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In Figures [5j [6j and [7] we show the partial decay widths as a function of the Higgs- 
boson mass for H — > eeqq, H — > euqq, and H — > qqqq, respectively. The upper plots 
show the predictions including both QCD and EW corrections. The lower plots depict 
the corrections relative to the lowest order. Besides the EW+QCD corrections, these 
plots include the EW and QCD corrections separately, the narrow-width approximation 
(NWA) and the improved Born approximation (IBA) as defined in Eqs. (7.5)-(7.7) and 
Eqs. (6.1)-(6.7), respectively, of Ref. [12]. We recall that the IBA for the partial decay 
widths includes leading effects such as corrections that are enhanced by factors G^m 2 
or G M M H , the Coulomb singularity for W pairs near their on-shell threshold, and the 
QCD correction to hadronically decaying gauge bosons. Apart from these effects, the 
IBA contains only one fitted constant for the WW- and ZZ-mediated channels each. Both 
in the WW- induced channel and in the ZZ-induced channel the EW corrections are very 
similar to the corresponding corrections for leptonic final states [12]. For moderate Higgs- 
boson mass, they are positive and below ~ 4% for decays via Z pairs. For the W-mediated 
decays the Coulomb singularity yields a large effect near the WW threshold and the EW 
corrections are in the range between 2% and 8% for moderate Higgs-boson mass. For all 
decays the EW corrections reach about 13% near Mh = 700 GeV. The thresholds for the 
on-shell decay of the Higgs boson into W bosons, Z bosons, and top quarks are manifest 
in the shape of the corrections. The QCD corrections amount to roughly a s /ir « 3.8% for 
semileptonic and 2a s /iT ~ 7.6% for hadronic final states and are practically independent of 
the Higgs-boson mass. For H — > eeqq the agreement between the full result and the NWA 
is (accidentally) at the per-mille level sufficiently above the ZZ threshold. For H — > evqq 
and H — > qqqq the NWA agrees with the full results within 1-2% above threshold. The 
IBA describes the full corrections within 2-3% for Mh < 400 GeV for all final states. 

In Section |3~T1 we defined a classification of QCD corrections for the four-quark final 
states. While only category (a), i.e. QCD corrections to gauge-boson decays, exists for 
the final states H — > qqQQ and H — > qq'QQ' (q ^ Q, Q'), all categories (a)-(d) contribute 
to H — > qqqq and H — > qqq'q' . Figure [S] shows the relative EW corrections and the subcon- 
tributions of the different categories of QCD corrections as a function of the Higgs-boson 
mass. The corrections to gauge-boson decays, i.e. category (a), make up practically all of 
the QCD part. Note that the contributions (b)-(d) are multiplied by a factor 10 in the 
plots. For Mh > 2Mw, these contributions are completely negligible. In this region they 
are suppressed by a factor (T v /M v ) 2 with respect to the leading contributions because 
they have two propagators less that can become resonant. Below the WW threshold this 
suppression becomes smaller but at Mh = 120 GeV the interference contribution is still 
rather small reaching only a few per mille. The largest corrections originate from interme- 
diate g*g* states [category (d)], because these corrections are proportional to a 2 a 2 rather 
than to a 3 a s as all other QCD corrections. 

4.3 Invariant-mass distributions 

In order to reconstruct the Higgs-decay events and in order to separate signal events 
from possible background events, distributions in the invariant mass of fermion pairs 
resulting from a W- or Z-boson decay should be investigated. On the l.h.s. of Figure [9] 
we show the invariant-mass distribution of the qq pair in the decay H — > eeqq including 
QCD and EW corrections for M H = 170 GeV and M H = 200 GeV, i.e. for one M H value 
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Figure 5: Partial decay width for H — > eeqq as a function of the Higgs-boson mass. The 
upper plots show the absolute prediction including QCD and EW corrections, and the 
lower plots show the relative size of the QCD and EW corrections separately, their sum 
(corrected) and the predictions of the NWA and the IBA. 
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Figure 6: Partial decay width for H — > euqq as a function of the Higgs-boson mass. The 
individual curves are defined as in Figure El 
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Figure 7: Partial decay width for H — > qqqq as a function of the Higgs-boson mass. The 
individual curves are defined as in Figure El 
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Figure 8: Comparison of the different QCD contributions defined in Section [3TT1 and the 
EW contribution to the corrections to the partial decay width for H — > qqqq as a function 
of the Higgs-boson mass. 



below and another above the on-shell threshold at 2Mz for Z-boson pairs. Above the 
threshold for the on-shell decay into a Z-boson pair there is a just a resonance around 
the Z-boson mass. Below the threshold only one Z boson can become resonant while the 
other Z boson is off shell. Hence, in addition to the peak around the Z-boson mass, the qq 
invariant-mass distribution shows an enhancement for M qq < Mh — M% « 80 GeV where 
the e + e~ pair can result from a resonant Z boson. 

The complete relative corrections to the distribution in the invariant mass of the qq 
pair and also of the e + e~ pair are shown on the r.h.s. of Figure [91 In addition, the 
QCD corrections to the M qq distribution are plotted separately; they are flat and amount 
to roughly 3.8%. Note that M qq actually is the total hadronic invariant mass resulting 
from the H — > eeqq decay, since we always recombine the qq(g) system to two jets. In a 
detailed experimental analysis a jet algorithm should be defined. Then, hard gluons can 
produce a separate jet and the QCD corrections need not be flat anymore. For such a 
study, the jet algorithm could simply be interfaced to our Monte Carlo program. The 
EW corrections reveal the same structure as discussed in the case of leptonic decays [12] 
shifting the peak position of the resonance. Close to the resonance and above, the EW 
corrections can reach 5—10%, below the resonance they become larger. The dominant 
effect is of photonic origin, leading to more pronounced corrections in the case of the 
leptonic invariant mass M e + e -, since the electric-charge factors are larger for leptons 
than for quarks. The way photons are treated has a strong impact on the corrections. 
By performing photon recombination, as defined in Section H~Tl we obtain collinear-safe 
observables. Thus, the corrections are of moderate size. However, for non-collinear-safe 
observables, i.e. if no photon recombination with leptons were performed, the corrections 
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Figure 9: Distribution in the invariant mass of the qq pair (l.h.s.) and relative EW+QCD 
corrections to the distributions in the invariant mass of the ee and qq pairs (r.h.s) in the 
decay H -> eeqq for M H = 170 GeV and M H = 200 GeV. For the distribution in M qq the 
relative QCD corrections are separately shown. 
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would be much larger because of mass- singular corrections proportional to aln(mi/Mn), 
as discussed in Ref. [12]. 

In Figure [TU] we show the distribution in the invariant mass of the qq pair and relative 
corrections to the distributions in the invariant mass of the ev and qq pairs in the decay 
H — ► evqq for Mh = 140 GeV and Mh = 170 GeV. Similarly to the decay H — > eegg there 
is a resonance around the W-boson mass and, for Mh < 2Mw, an additional enhancement 
for M qq < Mh — M\y « 60 GeV where the ev pair can become resonant. Also the relative 
corrections show the same characteristics. The corrections for the M eu distribution are 
somewhat smaller than for the M ee distribution in Figure O since the neutrino does not 
radiate photons. 

4.4 Angular distributions 

Angular distributions can be used to discriminate the Higgs-boson signal from the 
background or to study the properties of the Higgs boson. In Figure [TT] we show the 
distribution in the angle between the decay planes of the reconstructed Z bosons in the 
decay H — > eeqq in the rest frame of the Higgs boson. This angle can, for instance, 
be used to determine the parity of the Higgs boson [11]. Since the two jets cannot be 
distinguished, we show the distribution in the variable 

|cos^|= l( ,^ X ^f 1 kictlX 1 kict2 | )l , (4.2) 
|k had x ki||kjoti x kjotal 

which is symmetric with respect to the interchange of the jet momenta kj et i and kj Ct 2- 
Here the total hadronic momentum khad is equal to the sum of the two jet momenta, 
kj e ti + kj et2 , because we enforce 2-jet events, and ki is the momentum of the electron. For 
Mh = 200 GeV, both QCD and EW corrections are positive and about 4%. For Mh = 
170 GeV, the EW corrections are only about 2%. Both the EW and QCD corrections to 
this distribution are flat, in contrast to the EW corrections to the distribution in coscj)^ 
shown in Ref. [12] for analogous definitions of angles (p^ between the two planes defined 
by leptonically decaying Z bosons. This difference results from the fact that the sign of 
cos 0^ is only observable in the purely leptonic case. 

In the decay H — > evqq, angles between the electron and jets can be used for back- 
ground reduction [8]. In Figure [T2l we show the distribution in the angle between the 
electron and the W boson that is reconstructed from the qq pair in the rest frame of the 
Higgs boson and the corresponding relative QCD and EW corrections. The plot shows 
the well-known property that the electron is predominantly produced in the direction 
opposite to the hadronically decaying W boson. The QCD corrections are about 4% and 
the EW corrections at the level of 5%. The complete corrections can reach up to 12% 
depending on the value of the Higgs-boson mass. Since the EW corrections depend on 
the angle they distort the distribution by a few per cent. 
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Figure 10: Distribution in the invariant mass of the qq pair (l.h.s.) and relative EW+QCD 
corrections to the distributions in the invariant mass of the eu and qq pairs (r.h.s) in the 
decay H — > euqq for Mh = 140 GeV and Mh = 170 GeV. For the distribution in M qq the 
relative QCD corrections are separately shown. 
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Figure 11: Distribution in the angle between the Z — > ee and Z — > qq decay planes in the 
decay H — > eeqq (l.h.s.) and corresponding relative EW and QCD corrections (r.h.s.) for 
M H = 170 GeV and M H = 200 GeV. 



22 




Figure 12: Distribution in the angle between the electron and the W boson reconstructed 
from the qq pair (l.h.s.) and corresponding relative EW and QCD corrections (r.h.s.) in 
the decay H -> euqq for M H = 140 GeV and M H = 170 GeV. 
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5 Conclusions 



The decays of the Standard Model Higgs boson into four fermions via a W-boson or 
Z-boson pair lead to experimental signatures at the LHC and at a future e + e~ linear 
collider that are both important for the search for the Higgs boson and for studying 
its properties. In order to allow for adequate theoretical predictions for these decays, a 
Monte Carlo event generator is needed that properly accounts for the relevant radiative 
corrections. Prophecy4f is such an event generator which provides accurate predictions 
above, in the vicinity of, and below the WW and ZZ thresholds, owing to the use of the 
complex-mass scheme for the treatment of the gauge-boson resonances. 

While Prophecy4f originally contained only the electroweak corrections, in this 
paper we have included also the complete 0(a s ) QCD corrections. This allows to study 
precise predictions for all leptonic, semileptonic, and hadronic final states. 

The QCD corrections to the partial decay widths are dominated by the corrections to 
the gauge-boson decays and roughly given by a s /ir « 3.8% for semileptonic and 2a s /n ~ 
7.6% for hadronic final states. The electroweak corrections to the partial decay widths 
are very similar for leptonic, hadronic, and semileptonic final states. They are positive, 
typically amount to some per cent, increase with growing Higgs mass Mh, and reach 
about 8% at Mh ~ 500 GeV. In the on-shell (narrow-width) approximation for the 
intermediate gauge bosons, the correction is good within 1-2% of the partial widths for 
Higgs-boson masses sufficiently above the corresponding gauge-boson pair threshold, as 
long as the lowest-order prediction consistently includes the off-shell effects of the gauge 
bosons. For H — > WW — > 4/ the narrow-width approximation fails badly close to the WW 
threshold, because the instability of the W bosons significantly influences the Coulomb 
singularity near threshold. Only a calculation that keeps the full off-shellness of the 
W and Z bosons can describe the threshold regions properly. A simple improved Born 
approximation for the partial widths reproduces the full calculation within < 2-3% for 
Higgs-boson masses below 400 GeV. In this regime our complete calculation should have a 
theoretical uncertainty below 1%. For larger Higgs-boson masses we expect that unknown 
two-loop corrections that are enhanced by G M M^ deteriorate the accuracy. Finally, for 
Mft ^ 700 GeV it is well known that perturbative predictions become questionable in 
general. 

We have numerically investigated distributions for semileptonic final states where 
collinear photons are recombined and 2-jet events are forced. For angular and invariant- 
mass distributions the QCD corrections are flat and reflect the corresponding corrections 
to the integrated decay widths. For angular distributions, which can be used for back- 
ground reduction or the study of the quantum numbers of the Higgs boson, the electroweak 
corrections are of the order of 5-10% and, in general, distort the shapes. For invariant- 
mass distributions of fermion pairs, which are relevant for the reconstruction of the gauge 
bosons, well-known large photonic corrections show up and can exceed 10% depending on 
the treatment of photon radiation. 

This work completes the physics part of the Monte Carlo event generator Prophecy4f 
for H — > WW/ZZ — > 4/. It now includes the complete O(a) electroweak and 0(a s ) QCD 
corrections as well as corrections beyond 0(a) originating from heavy-Higgs effects and 
final-state photon radiation for all possible 4-fermion final states. Prophecy4f works 
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at the parton level and generates weighted events; unweighted event generation and an 
interface to parton showering will be addressed in the future. 
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